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Abstract —Considering a three-node multiple antenna relay 
system, this paper proposes a two-phase amplify-and-forward 
(AF) relaying protocol, which enables the autonomous relay to 
simultaneously harvest wireless power from the source informa¬ 
tion signal and from an energy signal conveyed by the destination. 
We first study this energy-flow-assisted (EFA) relaying in a single¬ 
input single-output (SISO) relay system and aim at maximizing 
the rate. By transforming the optimization problem into an 
equivalent convex form, a global optimum can be found. We 
then extend the protocol to a multiple antenna relay system. 
The relay processing matrix is optimized to maximize the rate. 
The optimization problem can be efficiently solved by eigenvalue 
decomposition, after linear algebra manipulation. It is observed 
that the benefits of the energy flow are interestingly shown only 
in the multiple antenna case, and it is revealed that the received 
information signal and the energy leakage at the relay can be 
nearly separated by making use of the signal space, such that 
the desired signal can be amplified with a larger coefficient. 

Index Terms —Energy harvesting, muftiple antenna reiay, reiay 
network, ampiify-and-forward (AF). 

I. Introduction 

During the uplink transmission in a sensor network, the 
relay may suffer from energy drain, while the destination, as 
a central processing unit, has access to a reliable source of 
energy. Motivated by this scenario, a joint wireless information 
and power transfer (JWIPT) 11 is investigated in an au¬ 
tonomous relay network, where the relay harvests energy from 
the incoming signal from the source to forward information to 
the destination but is also aided by a direct power transfer from 
the destination. The forms of wireless power harvesting and 
information relaying applied to current works on autonomous 
relays mainly focus on the power splitting (PS) relaying and 
the time switching (TS) relaying f2]|-|[5ll. Ref. |2j proposed 
a PS relaying (where the relay extracts power for forwarding 
from the source information signal) and a TS relaying (where 
the relay harvests power from an energy signal sent by the 
source and then relays source information in a time-division 
manner). Another TS relaying, where the energy signal is 
sent by the destination, was studied in (3). The PS relaying 
was also studied in multi-pair one-way relay networks E) 
and relay interference channels 0. In these works, the PS 
relaying reduces the received information power at the relay. 
The TS relaying consumes more timeslots, though the wireless 
power can be harvested in a dedicated timeslot. Thus, these 
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Fig. 1. JWIPT relay network. 

two methods may degrade the rate performance. To harvest 
sufficient power without consuming more timeslots, this paper 
proposes an energy-flow-assisted (EFA) two-phase amplify- 
and-forward (AF) one-way relaying protocol. 

In the proposed protocol, by exploiting PS, the relay can 
harvest power from both the source information signal and a 
concurrent energy flow (EF) at the destination (in the form 
of a WPT), as illustrated in Fig. |TJa). We first study the 
protocol for a single-input single-output (SISO) relay system 
and aim at maximizing the rate. By transforming the original 
nonconvex optimization problem into an equivalent convex 
form, a global optimum can be achieved. We then extend 
the protocol to a multiple antenna relay system. The relay 
processing matrix is optimized to maximize the rate. By linear 
algebra manipulation, the original problem is transformed 
into a generalized Rayleigh quotient problem which can be 
efficiently solved by eigenvalue decomposition. Simulation 
results reveal that only in the multiple antenna case, can the 
EFA relaying significantly outperform rate-wise the relaying 
without EF. It is shown that although the information signal 
to be forwarded is corrupted by the EF, the two signals can be 
nearly aligned with two orthogonal directions by manipulating 
the signal space, and the information signal can be amplified 
with a larger coefficient. 

The remainder of this paper is organized as follows. The 
system model is formulated in Section [II] Section [ill] then 
proposes the relaying protocols for the SISO and the multiple 
antenna relay systems. Section UlI-BI discusses the simulation 
results. Finally, conclusions are drawn in Section|V] Notations: 
Throughout the paper, matrices and vectors are in bold capital 
and bold lower cases, respectively. The notations (-) T , (•)*, 
(■) H , Tr{-}, vec(-), || • || p, and (•, •) represent the transpose, 
conjugate, conjugate transpose, trace, vectorization, Frobenius 
norm, Rronecker product, and the inner product, respectively. 

II. System Model 

Fig- Ea) shows the EFA two-phase relaying, where the 
single-antenna source S transmits information to the single¬ 
antenna destination D via a r-antenna wireless-powered half¬ 
duplex AF relay R, without the direct link between S and 
D. In phase 1, I) transmits an EF (i.e. an energy signal) in 
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the form of a power transfer to R, in order to enhance the 
amount of harvested power at the relay. Simultaneously, S 
transmits an information signal (i.e. an information flow) to R 
for forwarding. The relay R harvests power from both the EF 
and the information flow. In phase 2, the received signal at 
R is processed and forwarded thanks to the wireless power 
harvested in phase 1. The D-R, S-R, and R-D channels 
are respectively designated as h^D £ C rxl , £ C rxl , 

and \id,r £ C rxl , which are independent and identically 
distributed (i.i.d.) Rayleigh flat fading channels. Channel reci¬ 
procity is assumed such that h/j g = h/, D . Global channel 
state information (CSI) is supposed to be perfectly obtained at 
R. In phase 1, the relay exploits PS for simultaneous energy 
harvesting (EH) and information detecting (ID). Specifically, 
at each antenna of the relay, a fraction of the received power, 
denoted as the PS ratio p (where we consider the uniform 
PS), is conveyed to the EH receiver. The noise at the ID 
receiver at R and the receiver at D are respectively denoted 
by rig ~ CAf(0, er 2 I) and rip ~ CAf(0,a 2 ), while the noise 
at the EH receiver at R is small and neglected. 

In phase 1, the received signal at the EH receiver is given 
by yrj.EH = P 1/2 (^ir,dx d + h R>s x s ), where x D and x s 
respectively denote the EF from D and the information flow 
from S, and E{\xd\ 2 } = Pd and £{|a;s| 2 } = Ps, where Pn 
and Ps denote the power budgets at D and S. Assuming an 
RF-to-DC conversion efficiency of 1, the harvested wireless 
power at the EH receiver of R is given by 

Pr = £{||y,R,EH|| 2 } = p(\\^r.d\\ 2 Pd + Hh^slpPg). ( 1 ) 

Meanwhile, the received signal at the ID receiver is given by 
yti.iD = {i-p) 1/2 {h R ,sxs + h r^dXd) + n R . It shows that a 
part of the received EF (i.e. (1 — p) 1 // 2 h R:R >XD), referred to as 
an energy leakage, leaks into the ID receiver and contaminates 
the information signal. This leakage cannot be removed due 
to AF. In phase 2, y_R,iD is processed by the relay processing 
matrix F and forwarded from R: 

xr,id = F(1 ^ p) 1 ^ 2 (h R ,sXs + hjt’DXn) + Fn R . (2) 

The signal received at D is given by y' D = (1 — 

Py /2h -D,R Fh R,sXs+(l-p) 1/2 h^ R Fh R .DXD+h^ tR Fn R + 
no ■ Perfect knowledge of CSI and F at D is assumed. Thus, 
the self-interference (i.e. the term containing xd) is canceled, 
yielding 

y D = h^, fl F(l - p) 1/2 h R ,sxs + h^gFriR + no ■ (3) 

When Pd = 0 (i.e. xd = 0) as shown in Fig. [TJb), the 
forwarding power only comes from the information flow by 
PS. This relaying protocol is referred to as the relaying without 
EF. Note that although more power is available to the relay 
for forwarding in the EFA relaying, a part of the power used 
for forwarding is consumed to retransmit the energy leakage. 
In the following sections, we study the above two relaying 
protocols in both SISO and multiple antenna relay systems. 

III. Relaying Protocols 
A. The SISO System 

In the SISO case, each node of D, R , and S is equipped 
with one antenna, such that the D-R, S-R, and R-D channels 


are respectively designated as h RyR >, h R ^s, and /id, it- The 
relay processing matrix F also reduces to an amplification 
coefficient /. From 0, the achievable rate is obtained as 

Ci = 1/2 • log 2 (l + 7i), (4) 

where the coefficient 1/2 results from the half-duplex trans¬ 
mission, and the signal-to-noise ratio (SNR) 71 is given by 

^ _ (l-p)\h D , R \ 2 \f\ 2 \h R , s \ 2 Ps 

71 (1 + \h D>R \ 2 \f\ 2 )a 2 ’ 

where |/| 2 and p are variables to be optimized. Since the power 
of the amplified signal at R is no greater than the harvested 
wireless power, the power constraint at R is given by 

(l-p)Pi|/| 2 +^l /| 2 <PiP, ( 6 ) 

where pi = \h Rt s\ 2 Ps + !^-R,d| 2 Pd• It is found in (0 that 
C 1 increases as |/| 2 increases. When |/| 2 is maximized (i.e. 
Ci is maximized), equality holds in 0. Thus, 

\f\ 2 = PPi/((l~ P)Pl+&n)- I 7 ) 


The objective of the EFA relaying is maximizing the rate 
©. Since Ci is monotonically increasing, maximizing Ci is 
equivalent to maximizing 71 . Substituting © into 0. the 
optimization problem for the EFA relaying is formulated as 


_ P( 1 - P) _ 

p X (1 - p)pi + cr 2 + \h D ,R\ 2 Pip 
s.t. 0 < p < 1. 


( 8 a) 

(8b) 


Problem © is non-convex. However, the numerator of (f 8 ab . 
designated as <f>(p), is concave and nonnegative, while the 
denominator, designated as ijj(p), is convex and positive. Thus, 
© is a concave fractional programming. By introducing 
variables s = p/ip(p) and f = l/i/j(p), © can be equivalently 
transformed into a convex problem given by 


min —t<j>{s/t) (9a) 

s.t. ti/j(s/t) < 1, (9b) 

t > 0, 0 < s < t , (9c) 

which can be solved by the interior-point algorithm ©. Then, 
the optimal p* is obtained by p* = s*/t*, and the achievable 
rate © can be calculated with © and 0 . 

When Pd = 0, problem © boils down to the design 
problem of the relaying without EF. The optimal PS ratio can 
still be obtained by solving ©. 


B. The Multiple Antenna Relay System 

In the multiple antenna relay system (where r > 1), the 
system model is discussed in Section QI] According to © and 
©, the power of the forwarded signal and the SNR 72 of the 
received signal at D are respectively given by 

£{||x«, id || 2 } = (1 - p)Tv{F h F (h fl , s h % tS p s 

+ h R>D h* D P D )} +Tr{a 2 n F H F} (10) 

_ (l-p)P s |h^Fh fl , s | 2 
72 h£,*FF"h * DtR <jl+oV 
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Making use ofTr{ABC} = vec(A 7 r ) 7 r (C T (g)I)vec(B), ( ITOl ) 
can be written as 

^{||x fl ,i D || 2 }=vec(F) ff ((1 - p)Q/j(8)I+cr 2 l)vec(F),(12) 

where QK = h/i,S'h^ s Ps + hft rjiig jjPi). Similarly, thanks 
to vec(ABC) = (C T <g> A)vec(B), terms fl FhR g and 
h T DtR FF H h * in CD can be respectively rewritten as 

h D,fi Fh AS = vec ( h D.it Fh H,s) = (h R<s ® h D.fl) vec ( F ) 
and h^FF ff h^ fl = Tr{F ff h^h^F} = vec(F) ff (I 0 

^b.n^D ,it)vec(F). By defining f = vec(F), Qr = Q R 0 I, 

= [(hR i shR jg )®(h£) i R- 

h d , r )] T - an ^ J = = [I®(hz5,ith^ K )] 7 . {ED 

can be rewritten as 

72 = P s f H Kt/ (a 2 (f H Jf + 1)) , (13) 

where K = (1 — p) K. The rate of the multiple antenna relay 
system is obtained by C 2 = 1/2 ■ log 2 (l + 72 ). 

Similarly to the SISO system, maximizing the rate C 2 is 
equivalent to maximizing 72 . It is shown in the subsequent 
discussion that the coupled p and f are actually subject to an 
equality constraint (as shown in (115bb ). such that they cannot 
be optimized (or updated) alternatively. Therefore, as an initial 
study, in this SNR maximization problem, only f (i.e. the 
vectorization of the relay matrix) is optimized for simplicity. 
The best PS ratio p is exhaustively searched over all PS ratios 
to maximize the SNR in ©■ The power constraint at R is 
such that the transmit power (fl 2 l> is no greater than the power 
budget CD. Namely, 

(1 - p)f H Q R f + a 2 J H f < P R . (14) 

By decomposing f as f = f||f||, where f denotes the 
normalized direction vector of f, we find that the SNR 72 
in (fl~3l) increases as ||f || 2 increases. Hence, when (fT3l l is 
maximized, equality holds in otherwise, a larger ||f |J _ 2 can 
be found. Therefore, GS can be rewritten as ((1 — p)f H Q R i+ 
<r 2 f H f)/ P R = 1. With the above equation, the term i H Jf + 1 
in the denominator of (IT3l > can be rewritten as i H Jf, where 
J = J+((l— p)Qr+Ict 2 )/Pr. Hence, the SNR maximization 
problem can be formulated as 

max 7 ' = f ff Kf/ Jf) (15a) 

s.t. (1 - p)f H Q R f + alf H f = Pr , (15b) 

where (Q3j> is a generalized Rayleigh quotient which can 
then be reduced to a Rayleigh quotient by defining f = 
Ly x g = Ly 1 g||g||, where Lj is achieved by the Cholesky 
decomposition J = i// L j and g is the direction vector of g 
such that g H g = 1. Thus, problem ( 1 15 1 can be transformed 
into an equivalent form 


of (LJ^KL] 

solving (1 16bb with the achieved g* 
calculated with (THT i. 

To gain an insight into the benefits of multiple anten¬ 
nas, we rewrite (IT5ab such that f is reshaped as F. Given 
F = F/||F||r and its singular value decomposition F = 
USV 77 (where E = diag {\\,..., A r } for Ai > • • ■ > A r 
and Yh A 2 = 1 ), (1 1 5ab can be written as 7 ^ = (1 — 


ij . The optimum 


|g*|| can be obtained by 
Thereby, the rate can be 


E|hk B u 4 | 2 A? eiEIhgoV.pA 2 


P) I E K,R*i vf h_R,sAi| 2 /( 1 " + 


+ 


^E|h£ nV^A 2 


■ 2 ^-), where h p q denotes the normalized 
and v, denote the / th columns of U 
and V, respectively; d = (1 - p)Ps/{p\\h D , R \\ 2 ); e 2 = 

(1 - p)P D /(p\\h R , s \n P 2 = \\h R , D \\ 2 PD + ||hR, s || 2 Ps. 


P‘2 PP 2 

directions of h p>(? ; u/ 


Further, 7 ^ = (1 - p)\ J2 6 ,ie^ 1 ’* 6 .ie ie 2 ’*A J | 2 /( 


EC 


<uEG 2 


- + 




■||hR, S || 2_r 

P2 ■ _ P2 where £ M _e^ = (h* DR ,Ui), 

e302,i = (h = (h R ,D,Vi) such that 
Y^rn i = 1 f° r 171 = 1)2,3. The angle 9\p and the value 
denote the pseudo-angle and the cosine of the Hermitian 
angle (which is the angle between two complex vectors in 
the real vector space) between h* D R and u,, respectively 
Q. Likewise, £ 2 ,i> S 2i , £ 3 p, and 9 3 j denote corresponding 
pseudo-angles and cosines. Given fixed 9 m i , A, and £3 j, the 
value of 7 2 is increased, if the ordering of £1 j (and £ 2 , 1 ), 
i = 1 ,..., r, follows that of e.g. £i ? i > £l ,2 > • • • > £l,r 
for Ai > A 2 > • • • > A r . This implies that Ui and Vi should 
be close to h* D R and h R: s, respectively. Then, with fixed 
9m,i, A^ £i_i and £ 2i ,, to increase 7 2 , £ 3 ^ should be as large 
as possible, i.e. v r is close to h R ^. This v r can always be 
found by rotating v r around h^.s such that the angle between 
v r and does not change (i.e. £ 2 r is fixed). The above 

insight can also be verified by simulations. 

By exploiting maximum ratio combining (MRC) and maxi¬ 
mum ratio transmission (MRT), a simplified relay matrix can 
be designed as F' = yho ■ F ^ e baseline schemes with 

F' are referred to as EFA MRCMRT relaying and MRCMRT 
relaying without EF. The rate maximization problems (where 
// and p are the optimization variables) of these two MRCMRT 
protocols can be formulated as problems similar to ® and can 
also be equivalently transformed into convex forms similar to 
©. The details are omitted here, due to the space constraint. 
It is noteworthy that in the MRCMRT relaying without EF, 
F' = argmaxF 7 2 (p, F) ( 8 ], which indicates that the MR¬ 
CMRT structure is optimal for the scenario without EF. This 
is verified by Fig. [4] where the rate of the MRCMRT relaying 
without EF coincides with that of the relaying without EF 
(with an exhaustive search of p). 


IV. Simulation Results 


max g H (L/^KL/g (16a) 

i.llgll 

s.t. ||g|| 2 g ff [L 7 1 ] ff ((l- /9 )QR + ( T 2 l)L 7 1 g = PR.(16b) 

The optimal g* of problem (fl 6 l ) is equal to the dominant 
eigenvector of (L 7 1 ) KL 7 1 , such that the maximum value 

of the objective function (1 1 6ab is equal to the largest eigenvalue 


In the simulations, the channels are modeled as h p q = 
d pq h p q , where h M denotes the small-scale fading; d pq is 
the distance between nodes p and q, and d R s = 9 dr + d R s- 
In the simulations, dj^g — 10 m, cr 2 = 1 pW, and r = 4 in 
the multiple antenna relay system. The power budget settings 
consider the symmetric case of Pjj = 0.5 W and Pg = 0.1 W 
and the asymmetric case of P R = 5 W and Pg = 0.01 W. 
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(a) Rate as a function of power splitting ratio 
at a channel realization. 



(a) P D = 0.5 W, P s = 0.1 W. 



0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 
dnc/dnc ratio 


(a) P D = 0.5 W, P s = 0.1W. 



(b) Average rate as a function of the number 
of antennas. 

Fig. 2. The rate performance of the multiple 
antenna EFA relaying with dnn/dns = 0.5, 
P D = 0.5 W, and P s = 0.1 W. 




0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 
djdr, c ratio 


(b) P D = 5 W, P s = 0.01 W. (b) P D = 5 W, P s = 0.01 W. 

Fig. 3. Average rate as a function of Fig. 4. Average rate as a function of dpR/dps in 
dDR/dns i n the SISO relay system. the multiple antenna relay system. The best PS ratio is 

exhaustively searched among 0.01:0.01:0.99. 


Fig. |2(a)| shows that the rate curves of the EFA relaying 
and the relaying without EF are concave. This is because a 
low PS ratio results in a limited forwarding power budget, 
while a high PS ratio reduces the receive SNR at R. Both the 
above two factors can lead to a limited SNR at D. In Fig. 
|2(b)| the genie-aided EFA relaying refers to the EFA protocol 
where the energy leakage is ideally eliminated in the forwarded 
signal, i.e. the term containing xd in (0 is removed. Thus, no 
forwarding power is consumed for the energy leakage, and the 
rate of the genie-aided relaying is the upper bound, as shown in 
Fig. |2(b)| Since the relay beamforming of the EFA MRCMRT 
relaying does not address the energy leakage amplification 
problem, there is a gap between the rate of the EFA relaying 
and the EFA MRCMRT relaying. It is also observed that the 
rate of the EFA relaying can approach the ideal upper bound 
as the number of antennas increases. 

As shown in Fig. [3 the rate of the EFA relaying is negligibly 
higher than that of the relaying without EF in both the 
symmetric and the asymmetric cases. The reason lies in that 
the information signal and the energy leakage are amplified 
with the same amplification coefficient in the SISO relay, 
which means that the energy leakage may consume large 
amounts of the power for information forwarding. It is also 
observed that the rate of the EFA relaying nearly coincides 
with that of the relaying without EF when R is close to S 
in both the symmetric and asymmetric power budget cases. 
To figure out the reason, substituting 0 into 0 yields 71 = 

(1 - p)p\h DtR \' 2 \h Rj s\ 2 Ps/<^n/( l ~ P + <?l/Pl + \hD,R?P)- 
Thus, when R is close to S, may become very small 

such that the effect of Pjj on 71 is negligible. 

Rg. a illustrates that the rate of the EFA relaying can 
be much higher than that of the relaying without EF (and 
MRCMRT relaying without EF) in both the symmetric and 
asymmetric cases, since the forwarding power is enhanced 
by EF. The EFA relaying also outperforms rate-wise the EFA 
MRCMRT relaying. Different from the SISO system, the EFA 
relaying in multiple antenna system can utilize the signal 
space to project the main parts of the desired information 
and the energy leakage onto orthogonal directions such that 
the information signal to be forwarded can be amplified with 


a larger coefficient, while the EFA MRCMRT relaying does 
not suppress the amplification of the energy leakage. It is 
also observed that the gaps between the rate of the relaying 
protocols decrease as R moves towards the source. This is 
because the power of the EF received at R decreases. 

V. Conclusion 

In this paper, we have proposed an EFA two-phase relaying 
protocol. Protocols and optimization problems are designed 
and solved for both the SISO and the multiple antenna relay 
system. Simulation results show that EF is beneficial to the 
EFA relaying only in the presence of multiple antennas, 
because the information signal and the energy leakage can 
be nearly aligned with two orthogonal directions such that the 
desired signal can be amplified with a larger coefficient. The 
robust design addressing the problems of imperfect CSI and 
imperfect self-interference cancelation could be a future work. 
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